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Abstract 
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We give an overview of our precise calculation of the relic density of the lightest neutralino, in which we 
included relativistic Boltzmann averaging, subthreshold and resonant annihilations, and coannihilation processes 
with charginos and neutralinos. 



1. INTRODUCTION 

The lightest neutralino is one of the most 
promising candidates for the dark matter in the 
Universe. A linear combination of the superpart- 
ners of the neutral gauge and Higgs bosons, it is 
believed to be the lightest stable supersymmctric 
particle in the Minimal Supersymmetric exten- 
sion of the Standard Model (MSSM). 

In the near future, high precision measure- 
ments of the dark matter density may become 
possible from high resolution maps of the cosmic 
microwave background, and this may constrain 
supersymmetry. It is therefore of great interest 
to calculate the relic density of the lightest neu- 
tralino as accurately as possible. 

As a major step towards a complete and pre- 
cise calculation valid for all neutralino masses and 
compositions, we included for the first time Q all 
concomitant annihilations (coannihilations) be- 
tween neutralinos and charginos, properly treat- 
ing thermal averaging in presence of thresholds 
and resonances in the annihilation cross sections. 

2. FORMALISM 

Consider coannihilation of N supersymmetric 
particles with masses mi and statistical weights 
gi (first studied in ref. Normally, all heavy 
particles have time to decay into the lightest one, 
which we assume stable. Its final abundance is 
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then simply described by the sum of the densi- 
ties n = J2 n i- When the scattering rate of super- 
symmetric particles off the thermal background is 
much faster than their annihilation rate, n obeys 
the evolution equation || 



= -3Hn - (a cff v) (n 2 - n 2 cq ) 
with effective annihilation cross section 
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The numerator A is the total annihilation rate 
per unit volume at temperature T, and n oq is the 
total equilibrium density. Under the assumption 
of Boltzmann statistics (a good approximation for 
T < mi), we obtain [nj 
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Here Ki (x) is the modified Bessel function of the 
second kind of order i, and 
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where A(x, y, z) — x 2 + y 2 + z 2 — 2xy — 2xz — 2yz. 

W(s) is a Lorentz invariant annihilation rate 
per unit volume in which coannihilations appear 
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as thresholds at y/s equal to the sum of the masses 
of the coannihilating particles. The independence 
of W(s) on temperature is a remarkable calcula- 
tional advantage in presence of coannihilations: 
in fact it can be tabulated in advance, before tak- 
ing the thermal average and solving the density 
evolution equation. 

Eqs. (| jq) generalize the result of Gondolo and 
Gelmini [ff to coannihilations. 

3. RESULTS 

To explore a significant fraction of the MSSM 
parameter space |3| , we keep the number of theo- 
retical relations among the parameters to a min- 
imum. We assume GUT relations for gaug- 
ino masses, keep only the top and bottom tri- 
lincar soft supersymmetry-breaking parameters, 
and use a single mass parameter for the diago- 
nal entries in the sfermion mass matrices at the 
weak scale. We perform many different scans in 
parameter space, some general and some special- 
ized to interesting regions. We keep only models 
that satisfy the experimental constraints on the 
Z° width, on the b — > sj branching ratio, and on 
superpartner and Higgs boson masses. 
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Figure 1 . Neutralino relic density including neu- 
tralino and chargino coannihilations versus neu- 
tralino mass. The horizontal lines bound the cos- 
mologically interesting region 0.025 < il x h 2 < 1. 
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Figure 2. Ratio of the neutralino relic density 
with and without neutralino and chargino coan- 
nihilations versus neutralino composition. 



We obtain analytic expressions for the many 
Feynman diagrams contributing to the two-body 
cross sections at tree level for neutralino-neutral- 
ino, neutralino-chargino and chargino-chargino 
annihilation. Then for each set of model pa- 
rameters, we sum over particle polarizations and 
over initial and final states numerically and tab- 
ulate the annihilation rate W(s). Thermal av- 
eraging and integration of the density equation 
finally give the neutralino relic density Q, x h 2 = 
m x no/p C rit in units of the critical density p cx it- 

Fig. |l| shows the neutralino relic density fl x h 2 
versus the neutralino mass m x . Each point rep- 
resents a set of model parameters. It should be 
kept in mind that bands and holes in the point 
distributions are mere artifacts of our sampling 
in parameter space. 

The horizontal lines limit the cosmologically in- 
teresting region where the neutralino can consti- 
tute most of the dark matter in galaxies without 
violating the constraint on the age of the Uni- 
verse. We take it to be 0.025 < Q x h 2 < 1. 

The effect of coannihilations on the neutral- 
ino relic density is summarized in Fig. g which 
shows the ratio of the neutralino relic density with 
and without coannihilations versus the neutralino 
composition Z g /(l — Z g ) (here Z g is the gaugino 
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fraction) . 

Coannihilation processes are important not 
only for light higgsino-like neutralinos, as pointed 
out before in approximate calculations ||, but 
also for heavy higgsinos and for mixed and 
gaugino-like neutralinos. Indeed, coannihilations 
should be included whenever < 2|Mi|, inde- 
pendently of the neutralino composition. When 
|^| ~ | Mi |, coannihilations can increase or de- 
crease the relic density in and out of the cosmo- 
logically interesting region. 

Fig. |^ shows the cosmologically interesting re- 
gion in the neutralino mass-composition plane, 
before and after including coannihilations. This 
region is limited to the left by accelerator con- 
straints, to the right by Vl x h 2 < 1, and below 
and above by incompleteness in our survey of pa- 
rameter space, except for the hole at 85-450 GeV 
where VL x h 2 < 0.025. 

The main effect of coannihilations is to shift the 
higgsino region to higher masses. In particular 
the cosmological upper bound on the neutralino 
mass changes from 3 to 7 TeV. Differently from 
previous approximate results ||] , there is a cosmo- 
logically interesting window of light higgsino-like 
neutralinos with masses around 75 GeV. 

We conclude that if coannihilations are prop- 
erly included, the neutralino is a good dark mat- 
ter candidate whether it is light or heavy, and 
whether it is higgsino-like, mixed or gaugino-like. 

P. Gondolo thanks P. Salati for bringing ref. || 
to his attention. 
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Figure 3. Neutralino masses m x and composi- 
tions Z g /(1 — Zg) for cosmologically interesting 
models (a) with and (b) without coannihilations. 
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